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Abstract 

This paper provides a survey of known results and open problems for the 
two-type Richardson model, which is a stochastic model for competition 
on In its simplest formulation, the Richardson model describes the 
evolution of a single infectious entity on but more recently the dy¬ 
namics have been extended to comprise two competing growing entities. 
For this version of the model, the main question is whether there is a 
positive probability for both entities to simultaneously grow to occupy 
infinite parts of the lattice, the conjecture being that the answer is yes if 
and only if the entities have the same intensity. In this paper attention 
focuses on the two-type model, but the most important results for the 
one-type version are also described. 
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1 Introduction 

Consider an interacting particle system in which, at any time t, each site 
X G is in either of two states, denoted by 0 and 1. A site in state 0 flips to 
a 1 at rate proportional to the number of nearest neighbors in state 1, while 
a site in state 1 remains a 1 forever. We may think of sites in state 1 as being 
occupied by some kind of infectious entity, and the model then describes the 
propagation of an infection where each infected site tries to infect each of its 
nearest neighbors on Z'^ at some constant rate A > 0. More precisely, if at 
time t a vertex x is infected and a neighboring vertex y is uninfected, then, 
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conditional on the dynamics up to time t, the probability that x infects y 
during a short time window (t, t + h) is \h + o{h). Here and in what follows, 
sites in state 0 and 1 are referred to as uninfected and infected respectively. 
This is the intuitive description of the model; a formal definition is given in 
Section [2l 

The model is a special case of a class of models introduced by Richardson 
(1973), and is commonly referred to as the Richardson model. It has several 
cousins among processes from mathematical biology, see e.g. Eden (1961), 
Williams and Bjerknes (1972) and Bramson and Griffeath (1981). The model 
is also a special case of so called first-passage percolation, which was introduced 
in Hammersley and Welsh (1965) as a model for describing the passage of a 
fluid through a porous medium. In first-passage percolation, each edge of the 
Z'^-lattice is equipped with a random variable representing the time it takes 
for the fluid to traverse the edge, and the Richardson model is obtained by 
letting these passage times be i.i.d. exponential. 

Since an infected site stays infected forever, the set of infected sites in the 
Richardson model increases to cover all of as t —>■ oo, and attention focuses 
on how this set grows. The main result is roughly that the infection grows 
linearly in time in each fixed direction and that, scaled by a factor 1/t, the 
set of infected points converges to a non-random asymptotic shape as t —>■ oo. 
To prove that the growth is linear in a fixed direction involves Kingman’s 
subadditive ergodic theorem - in fact, the study of first-passage percolation 
was one of the main motivations for the development of subadditive ergodic 
theory. That the linear growth is preserved when all directions are considered 
simultaneously is stated in the celebrated shape theorem (Theorem 12.11 in 
Section [3]) which originates from Richardson (1973). 

Now consider the following extension of the Richardson model, known as 
the two-type Richardson model and introduced in Haggstrom and Pemantle 
(1998). Instead of two possible states for the sites there are three states, 
which we denote by 0, 1 and 2. The process then evolves in such a way that, 
for i = 1,2, a site in state 0 flips to state i at rate A* times the number of 
nearest neighbors in state i and once in state 1 or 2, a site remains in that 
state forever. Interpreting states 1 and 2 as two different types of infection 
and state 0 as absence of infection, this gives rise to a model describing the 
simultaneous spread of two infections on Z'’*. To rigorously define the model 
requires a bit more work; see Section [3l In what follows we will always assume 
that d> 2] the model makes sense also for d = 1 but the questions considered 
here become trivial. 

A number of similar extensions of (one-type) growth models to (two-type) 
competition models appear in the literature; see for instance Neuhauser (1992), 
Durrett and Neuhauser (1997), Kordzakhia and Talley (2005) and Ferrari et al. 
(2006). These tend to require somewhat different techniques, and results tend 
not to be easily translated from these other models to the two-type Richardson 
model (and vice versa). Closer to the latter are (non-Markovian) competition 
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models based on first-passage percolation models with non-exponential pas¬ 
sage time variables - Caret and Marchand (2005), Hoffman (2005:1), Hoffman 
(2005:2), Caret and Marchand (2006), Couere (2007), Pimentel (2007) - and 
a certain continuum model - Deijfen et al. (2004), Deijfen and Haggstrom 
(2004), Couere (2007). For ease of exposition, we shall not consider these 
variations even in cases where results generalize. 

The behavior of the two-type Richardson model depends on the initial 
configuration of the infection and on the ratio between the intensities Ai and 
A 2 of the infection types. Assume first, for simplicity, that the model is started 
at time 0 from two single sites, the origin being type 1 infected and the site 
(1,0,..., 0) next to the origin being type 2 infected. Three different scenarios 
for the development of the infection are conceivable: 

(a) The type 1 infection at some point completely surrounds type 2, thereby 
preventing type 2 from growing any further. 

(b) Type 2 similarly strangles type 1. 

(c) Both infections grow to occupy infinitely many sites. 

It is not hard to see that, regardless of the intensities of the infections, out¬ 
comes (a) and (b) where one of the infection types at some point encloses the 
other have positive probability regardless of Ai and A 2 . This is because each 
of (a) and (b) can be guaranteed through some finite initial sequence of infec¬ 
tions. In contrast, scenario (c) - referred to as infinite coexistence - can never 
be guaranteed from any finite sequence of infections, and is therefore harder 
to deal with: the main challenge is to decide whether, for given values of the 
parameters Ai and A 2 , this event (c) has positive probability or not. Intu¬ 
itively, infinite coexistence represents some kind of power balance between the 
infections, and it seems reasonable to suspect that such a balance is possible 
if and only if the infections are equally powerful, that is, when Ai = A 2 . This 
is Conjecture 13.31 in Section [Sj which goes back to Haggstrom and Pemantle 
(1998), and, although a lot of progress have been made, it is not yet fully 
proved. We describe the state of the art in Sections 0] and [5l 

As mentioned above, apart from the intensities, the development of the in¬ 
fections in the two-type model also depends on the initial state of the model. 
However, if we are only interested in deciding whether the event of infinite 
coexistence has positive probability or not, it turns out that, as long as the 
initial configuration is bounded and one of the sets does not completely sur¬ 
round the other, the precise configuration does not matter, that is, whether 
infinite coexistence is possible or not is determined only by the relation be¬ 
tween the intensities. This is proved in Deijfen and Haggstrom (2006:1); see 
Theorem 13.21 in Section [3] for a precise formulation. Of course one may also 
consider unbounded initial configurations. Starting with both infection types 
occupying infinitely many sites means - apart from in very labored cases - 
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that they will both infect infinitely many sites. A more interesting case is 
when one of the infection types starts from an inhnite set and the other one 
from a finite set. We may then ask if outcomes where the finite type infects 
infinitely many sites have positive probability or not. This question is dealt 
with in Deijfen and Haggstrom (2007), and we describe the results in Section 

E 

The dynamics of the two-type Richardson model is deceptively simple, and 
yet gives rise to intriguing phenomena on a global scale. In this lies a large 
part of the pleasure indicated in the title. Furthermore, proofs tend to involve 
elegant probabilistic techniques such as coupling, subadditivity and stochastic 
comparisons, adding more pleasure. The pain alluded to (which by the way 
is not so severe that it should dissuade readers from entering this field) comes 
from the stubborn resistance that some of the central problems have so far put 
up against attempts to solve them. A case in point is the “only if” direction 
of the aforementioned Conjecture 13.31 saying that infinite coexistence starting 
from a bounded initial configuration does not occur when Ai 7 ^ A 2 . 

2 The one-type model 

As mentioned in the introduction, the one-type Richardson model is equivalent 
to first-passage percolation with i.i.d. exponential passage times. To make the 
construction of the model more precise, first define E^d as the edge set for 
the lattice (i.e., each pair of vertices x,y € Z'^ at Euclidean distance 1 
from each other have an edge e € connecting them). Then attach i.i.d. 
non-negative random variables {T{e)}e£E^a edges. We take each r(e) 

to be exponentially distributed with parameter A > 0, meaning that 

P(r(e) > t) = exp(—Af) 
for all t > 0. For x, y € Z'^, define 

r(x,y) = inf^r(e) ( 1 ) 

eer 

where the infimum is over all paths F from x to y. The Richardson model 
with a given set Sq C Z'^ of initially infected sites is now defined by taking the 
set St of sites infected at time t to be 

St = {x G : T{y, x) <t for some y G 5o}. (2) 

It turns out that the infimum in ([1]) is a.s. a minimum and attained by a unique 
path. That St grows in the way described in the introduction is a consequence 
of the memoryless property of the exponential distribution: for any s,t > 0 
we have that P(T(e) > s + t |T(e) > s) = exp(—Af). 

Note that for any x,y,z G Z'^ we have T{x,y) < T{x,z) +T{z,y). This 
subadditivity property opens up for the use of subadditive ergodic theory in 
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analyzing the model. To formulate the basic result, let T{x) be the time when 
the point x € Z'^ is infected when starting from a single infected site at the 
origin and write n = (n, 0,..., 0). It then follows from the subadditive ergodic 
theorem - see e.g. Kingman (1968) - that there is a constant such that 
T{n)/n —>■ fix almost surely and in Li as n —> oo. Furthermore, a simple time 
scaling argument implies that fix = Xfii and hence, writing fii = fi, we have 
that 


lim -= Xfi a.s. and in Li. 

n—>-oo n 

The constant fi indicates the inverse asymptotic speed of the growth along the 
axes in a unit rate process and is commonly referred to as the time constant. 
It turns out that fi > 0, so that indeed the growth is linear in time. Similarly, 
an analog of ([3]) holds in any direction, that is, for any x € Z'^, there is a 
constant fi{x) > 0 such that T{nx)/n —>■ Xfi{x). The infection hence grows 
linearly in time in each hxed direction and the asymptotic speed of the growth 
in a given direction is an almost sure constant. 

We now turn to the shape theorem, which asserts roughly that the linear 
growth of the infection is preserved also when all directions are considered 
simultaneously. More precisely, when scaled down by a factor 1/t the set St 
converges to a non-random shape A. To formalize this, let St C be a 
continuum version of St obtained by replacing each x G by a unit cube 
centered at x. 


Theorem 2.1 (Shape Theorem) There is a compact convex set A such 
that, for any e > 0, almost surely 

(1 - e)XA c y C (1 + e)XA 

for large t. 

In the above form, the shape theorem was proved in Kesten (1973) as an 
improvement on the original “in probability” version, which appears already 
in Richardson (1973). See also Cox and Durrett (1988) and Boivin (1991) 
for generalizations to first-passage percolation processes with more general 
passage times. Results concerning fluctuations around the asymptotic shape 
can be found, e.g., in Kesten (1993), Alexander (1993) and Newman and Piza 
(1995), and, for certain other passage time distributions, in Benjamini et al. 
(2003). 

Working out exactly, or even approximately, what the asymptotic shape 
A is has turned out to be difficult. Obviously the asymptotic shape inherits 
all symmetries of the Z'^ lattice - invarince under reflection and permutation 
of coordiante hyperplanes - and it is known to be compact and convex, but, 
apart from this, not much is known about its qualitative features. These dif¬ 
ficulties with characterizing the shape revolve around the fact that Z'^ is not 
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rotationally invariant, which causes the growth to behave differently in differ¬ 
ent directions. For instance, simulations on W? indicate that the asymptotic 
growth is slightly faster along the axes as compared to the diagonals. There 
is however no formal proof of this. 

Before proceeding with the two-type model, we mention some work con¬ 
cerning properties of the time-minimizing paths in ([1]), also known as geodesics. 
Starting at time 0 with a single infection at the origin 0, we denote by r(x) 
the (unique) path F for which the infimum r(0,x) in © is attained. Define 
'h = U^^^dr(x), making T a graph specifying which paths the infection ac¬ 
tually takes. It is not hard to see that T is a tree spanning all of and 
hence there must be at least one semi-infinite self-avoiding path from the ori¬ 
gin (called an end) in T. The issue of whether T has more than one end was 
noted by Haggstrom and Pemantle (1998) to be closely related to the issue of 
infinite coexistence in the two-type Richardson model with Ai = A 2 : such infi¬ 
nite coexistence happens with positive probability starting from a finite initial 
configuration if and only if T has at least two ends with positive probability. 

We say that an infinite path xi,X 2 , ■ ■ ■ has asymptotic direction x if Xk/\xk \ —>■ 
X as A: —7> 00 . In d = 2, it has been conjectured that every end in T has an 
asymptotic direction and that, for every x € M^, there is at least one end (but 
never more than two) in T with asymptotic direction x. In particular, this 
would mean that T has uncountably many ends. For results supporting this 
conjecture, see Newman (1995) and Newman and Licea (1996). In the former 
of these papers, the conjecture is shown to be true provided an unproven but 
highly plausible assumption on the asymptotic shape A, saying roughly that 
the boundary is sufficiently smooth. See also Lalley (2003) for related work. 

Results not involving unproven assumptions are comparatively weak: The 
coexistence result of Haggstrom and Pemantle (1998) shows for d = 2 that 'k 
has at least two ends with positive probability. This was later improved to T 
having almost surely at least 2d ends, by Hoffman (2005:2) for d = 2 and by 
GouCTe (2007) for higher dimensions. 

3 Introducing two types 

The definition of the two-type Richardson model turns out to be simplest in the 
symmetric case Ai = A 2 , where the same passage time variables {r(e)}ee£;^^ 
as in the one-type model can be used, with A = Ai = A 2 . Suppose we start 
with an initial configuration (5g, Sq) of infected sites, where Sq C Z'^ are those 
initially containing type 1 infection, and ^g C Z'^ are those initially containing 
type 2 infection. We wish to define the sets S} and of type 1 and type 2 
infected sites for all t > 0. To this end, set Sq = S'g U 5g, and take the 
set St = St U St of infected sites at time t to be given by precisely the same 
formula ([21) as in the one-type model; a vertex x G is then assigned infection 
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1 or 2 depending on whether the y £ Sq for which 

inf{r(y,x) : 2/ G S’o} 

is attained is in Sq or Sq. 

As in the one-type model, it is a straightforward exercise involving the 
memoryless property of the exponential distribution to verify that (S^, S^)t>o 
behaves in terms of infection intensities as described in the introduction. 

This construction demonstrates an intimate link between the one-type and 
the symmetric two-type Richardson model: if we watch the two-type model 
wearing a pair of of glasses preventing us from distinguishing the two types 
of infection, what we see behaves exactly as the one-type model. The link 
between infinite coexistence in the two-type model and the number of ends in 
the tree of infection 'k of the one-type model claimed in the previous section 
is also a consequence of the construction. 

In the asymmetric case Ai ^ A 2 , the two-type model is somewhat less 
trivial to define due to the fact that the time it takes for infection to spread 
along a path depends on the type of infection. There are various ways to deal 
with this, one being to assign, independently to each e G E^d, two independent 
random variables ri(e) and r 2 (e), exponentially distributed with respective 
parameters Ai and A 2 , representing the time it takes for infections 1 resp. 2 to 
traverse e. Starting from an intial configuration (S'q,5q), we may picture the 
infections as spreading along the edges, taking time ri(e) or r 2 (e) to cross e 
depending on the type of infection, with the extra condition that once a vertex 
becomes hit by one type of infection it becomes inaccessible for the other type. 
This is intuitively clear, but readers with a taste for detail may require a more 
rigorous definition, which however we refrain from here; see Haggstrom and 
Pemantle (2000) and Deijfen and Haggstrom (2006:1). 

We now move on to describing conjectures and results. Write Gi for the 
event that type i infects infinitely many sites on and define G = Gi fl G 2 . 
The question at issue is: 

Does G have positive probability? (4) 

A priori, the answer to this question may depend both on the initial configu¬ 
ration - that is, on the choice of the sets Sq and 5 q - and on the ratio between 
the infection intensities Ai and A 2 . However, it turns out that, if we are not 
interested in the actual value of the probability of G, but only in whether it 
is positive or not, then the initial configuration is basically irrelevant, as long 
as neither of the initial sets completely surrounds the other. This motivates 
the following definition. 

Definition 3.1 Let .^1 and ^2 be two disjoint finite subsets of . We say 
that one of the sets (f,i) strangles the other (f^j) if there exists no infinite self¬ 
avoiding path in that starts at a vertex in f,j and that does not interseet f,i. 
The pair (^ 1 ,^ 2 ) is said to be fertile if neither of the sets strangles the other. 
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Now write for the distribution of a two-type process started from 

Sq = and Sq = ^ 2 - We then have the following result. 

Theorem 3.2 Let (^ 1 ,^ 2 ) and (^^,^ 2 ) fertile pairs of disjoint finite 

subsets ofU^, where d>2. For all choices o/(Ai,A 2 ), we have 

>0 4^ ^ 

For connected initial sets and ^2 and d = 2, this result is proved in 
Haggstrom and Pemantle (1998). The idea of the proof in that case is that, 
by controlling the passage times of only finitely many edges, two processes 
started from (^ 1 ,^ 2 ) and (^ 1 )^ 2 ) respectively can be made to evolve to the 
same total infected set after some finite time, with the same configuration of 
the infection types on the boundary. Coupling the processes from this time 
on and observing that the development of the infections depends only on the 
boundary configuration yields the result. This argument however breaks down 
when the initial sets are not connected (since it is then not sure that the same 
boundary configuration can be obtained in the two processes) and it is unclear 
whether it applies for d > 3. Theorem l3.2l is proved in full generality in Deijfen 
and Haggstrom (2006:1), using a more involved coupling construction. 

It follows from Theorem 13.21 that the answer to (j4]) depends only on the 
value of the intensities Ai and A 2 . Hence it is sufficient to consider a process 
started from S'q = 0 and 5q = 1 (recall that n = (n, 0,... ,0)), and in this 
case we drop subscripts and write P^i ’^2 fQj, Also, by time-scaling, we 

may assume that Ai = 1. The following conjecture, where we write A 2 = A, 
goes back to Haggstrom and Pemantle (1998). 

Conjecture 3.3 In any dimension d > 2, we have that P^’^{G) > 0 if and 
only if X = 1. 

The conjecture is no doubt true, although proving it has turned out to be a 
difficult task. In fact, the “only if” direction is not yet fully established. In 
the following two sections we describe the existing results for A = 1 and X 1 
respectively. 

4 The case A = 1 

When A = 1, we are dealing with two equally powerful infections and Conjec¬ 
ture 13.31 predicts a positive probability for infinite coexistence. This part of 
the conjecture has been proved: 

Theorem 4.1 // A = 1, we have, for any d>2, that P^’^(G) > 0. 


This was first proved in the special case d = 2 by Haggstrom and Pemantle 
(1998). That proof has a very ad hoc flavor, and heavily exploits not only 
the two-dimensionality but also other specific properties of the square lattice, 
including a lower bound on the time constant /r in ([3|) that just happens 
to be good enough. When eventually the result was generalized to higher 
dimensions, which was done simultaneously and independently by Caret and 
Marchand (2005) and Hoffman (2005:1), much more appealing proofs were 
obtained. Yet another distinct proof of Theorem 14.11 was given by Deijfen 
and Haggstrom (2007). All four proofs are different, though if you inspect 
them for a smallest common denominator you find that they all make critical 
use of the fact that the time constant /i is strictly positive. We will give the 
Garet-Marchand proof below. In Hoffman’s proof ergodic theory is applied to 
the tree of infection ik and a so-called Busemann function which is shown to 
exhibit contradictory behavior under the assumption that infinite coexistence 
has probability zero. The Deijfen-Haggstrom proof proceeds via the two-type 
Richardson model with certain infinite initial configurations (cf. Section [6]) . 

Proof of Theorem I4.lt The following argument is due to Caret and Marc¬ 
hand (2005), though our presentation follows more closely the proof of an 
analogous result in a continuum setting in Deijfen and Haggstrom (2004) - a 
paper that, despite the publication dates, was preceded by and also heavily 
influenced by Caret and Marchand (2005). 

Fix a small e > 0. By Theorem 13.21 we are free to choose any finite 
starting configuration we want, and here it turns out convenient to begin with 
a single type 1 infection at the origin 0, and a single type 2 infection at a 
vertex n = (n, 0,..., 0), where n is large enough so that 

(i) i?[r(0,n)] < (l-|-e)n/r, and 

(ii) P(r(0, n) < (1 — e)n/r) < e; 

note that both (i) and (ii) hold for n large enough due to the asymptotic 
speed result ([3]). The reader may easily check, for later reference, that (i) and 
(ii) together with the nonnegativity of T(0,n) imply for any event B with 
P{B) = a that 


Next comes an important telescoping idea: for any positive integer k we have 

.E[T(0,/cn)] = .F[r(0,n)]+.E[r(0,2n)-r(0,n)]+.E[r(0,3n)-T(0,2n)] 

+ ... + kn) - r(0, {k - l)n)] . 

Since limfc_).oo A:“^F1[T(0, A;n)] = n/i, there must exist arbitrarily large k such 
that 

E[T{0, {k + l)n) — T(0, A;n)] > (1 — e)nfj,. 


H[r(0,n)hH] 


3e 


1 — a 


n/r. 


(5) 
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By taking m = /cn, and by translation and reflection invariance, we may 
deduce that 

-E[T(n, —m) — T(0, —m)] > (1 — e)n/r (6) 

for some arbitrarily large m. We will pick such an m; how large will soon be 
specified. 

The goal is to show that P{G) > 0, so we may assume for contradiction 
that P{G) = 0. By symmetry of the initial configuration, we then have that 
P{Gi) = P{G 2 ) = This implies that 

lim P(—m gets infected by type 2) = lim P(T(n, —m) < r(0, —m)) = — 

m—>-co m^oo 2 

so let us pick m in such a way that 

P(r(n,-m) <T(0,-m)) > i (7) 

while also ([6]) holds. Write B for the event in ([7]). The expectation E[T{n, —m) — 
T(0, —m)] may be decomposed as 


£’[r(n, —m) — T(0, —m)] = 


< 

< 

< 

< 


E[T{n, -m) - T(0, -m) \B]P{B) 
+E[T{n, -m) - r(0, -m) \^B]P{^B) 
E[T{n, -m) - T(0, -m) \^B]P{^B) 

^P[r(n,-m) - r(0,-m) hP] 
^P[r(n,0)hP] 

3 . 

-(1 + 4e)n/r 


where the second-to-last inequality is due to the triangle inequality T(n, —m) < 
T(n,0) +T(0,—m), and the last one uses (l5|). For small e, this contradicts 
(| 6 |l, so the proof is complete. □ 


5 The case A 7 ^ 1 

Let us move on to the case when A 7 ^ 1, that is, when the type 2 infection 
has a different intensity than type 1. It then seems unlikely that the kind of 
equilibrium which is necessary for infinite coexistence to occur would persist 
in the long run. However, this part of Conjecture 13.31 is not proved. The best 
result to date is the following theorem from Haggstrom and Pemantle (2000). 

Theorem 5.1 For any d>2, we have P^’^{G) = 0 for all but at most count¬ 
ably many values of X. 
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We leave it to the reader to decide whether this is a very strong or a very 
weak result: it is very strong in the sense of showing that infinite coexistence 
has probability 0 for (Lebesgue)-almost all A, but very weak in the sense that 
infinite coexistence is not ruled out for any given A. 

The result may seem a bit peculiar at first sight and we will spend some 
time explaining where it comes from and where the difficulties arise when one 
tries to strengthen it. Indeed, as formulated in Conjecture l.'f.,'Il the belief is 
that the set {A : > 0} in fact consists of the single point A = 1, but 

Theorem 15. II only asserts that the set is countable. 

First note that, by time-scaling and symmetry, we have 
and hence it is enough to consider A < 1. An essential ingredient in the proof 
of Theorem IQ is a coupling of the two-type processes obtained 

by associating two independent exponential mean 1 variables ri(e) and T^ie) 
to each edge e € and then letting the type 2 passage time at parameter 
value A be given by T 2 (e) = A“^T 2 (e) and the type 1 time (for any A) by Ti(e). 
Write Q for the probability measure underlying this coupling and let G^ be 
the event that infinite coexistence occurs at parameter value A. Theorem 15.11 
is obtained by showing that 

with Q-probability 1 the event G^ occurs , . 

for at most one value of A G (0,1]. ^ 

Hence, Q{G^) can be positive for at most countably many A, and Theorem 
15.11 then follows by noting that P^’^{G) = Q{G^). 

But why is ([8]) true? Let Gf be the event that the type i infection grows 
unboundedly at parameter value A. Then the coupling defining Q can be shown 
to be monotone in the sense that G^ is decreasing in A - that is, if G^ occurs 
then Gi occurs for all A^ < A as well - and G 2 is increasing in A. This kind 
of monotonicity of the coupling is crucial for proving ([8]), as is the following 
result, which asserts that, on the event that the type 2 infection survives, the 
total infected set in a two-type process with distribution P^’^, where A < 1, 
grows to a first approximation like a one-type process with intensity A. More 
precisely, the speed of the growth in the two-type process is determined by the 
weaker type 2 infection type. We take SI to denote the union of all unit cubes 
centered at points in SI and A is the limiting shape for a one-type process 
with rate 1. 

Theorem 5.2 Consider a two-type process with distribution gome 

A < 1. On the event G 2 we have, for any e > 0, that almost surely 

cl I I c2 

(1 - £)XA c — - C (1 -h e)AA 

for large t. 
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Theorem [Q follows readily from this result and the monotonicity properties 
of the coupling Q. Indeed, fix e > 0 and suppose occurs. Then Theorem 
15.21 guarantees that on level A the type 1 infection is eventually contained in 
(1 + e)XtA, a conclusion that extends to all X' > A, because increasing the 
type 2 infection rate does not help type 1. On the other hand, for any X' > X 
we get on level A' that the union of the two infections will - again by Theorem 
15.21 - eventually contain (1 — e)X'tA, so by taking e sufficiently small we see 
that the type 1 infection is strangled on level A', implying ([8|), and Theorem 
15.11 follows. 

We will not prove Theorem l5.2l but mention that the hard work in proving 
it lies in establishing a certain key result (Proposition 2.2 in Haggstrom and 
Pemantle (2000)) that asserts that if the strong infection type reaches outside 
(1 + e)XtA infinitely often, then the weak type is doomed. The proof of this 
uses geometrical arguments, the most important ingredient being a certain 
spiral construction, emanating from the part of the strong of infection reach¬ 
ing beyond (1 -|- e)XtA, and designed to allow the strong type to completely 
surround the weak type before the weak type catches up from inside. 

How would one go about to strengthen Theorem 15.11 and rule out inhnite 
coexistence for all A 1? One possibility would be to try to derive a contra¬ 
diction with Theorem 15.21 from the assumption that the strong infection type 
grows unboundedly. For instance, intuitively it seems likely that the strong 
type occupying a positive fraction of the boundary of the infected set would 
cause the speed of the growth to exceed the speed prescribed by the weak 
infection type. This type of argument is indeed used in Caret and Marchand 
(2007) to show, for d = 2, that on the event of infinite coexistence the fraction 
of infected sites occupied by the strong infection will tend to 0 as t —>■ oo. This 
feels like a strong indication that infinite coexistence does not happen. 

Another approach to strengthening Theorem 15.11 in order to obtain the 
only-if direction of Conjecture 13.31 is based on the observation that, since co¬ 
existence represents a power balance between the infections, it is reasonable 
to expect that decreases as A moves away from 1. We may formulate 

that intuition as a conjecture: 

Conjecture 5.3 For the two-type Richardson model on TA with d > 2, we 
have, for A < A' G (0,1], that P^’^{G) < pi’^'(G). 

A confirmation of this conjecture would, in combination with Theorem 15.11 
clearly establish the only-if direction of Conjecture 13.31 If P^’^(G) > 0 for 
some A < 1, then, according to Conjecture 15.31 we would have (G) > 0 
for all X' G (A, 1] as well. But the interval (A, 1] is uncountable, yielding a 
contradiction to Theorem EH 

Although Conjecture 15.31 might seem close to obvious, it has turned out to 
be very difficult to prove. A natural first attempt would be to use coupling. 
Consider for instance the coupling Q described above. As pointed out, the 
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events and G 2 that the individual infections grow unboundedly at param¬ 
eter value A are then monotone in A, but one of them is increasing and the 
other is decreasing, so monotonicity of their intersection does not follow. 
Hence more sophisticated arguments are needed. 

Observing how our colleagues react during seminars and corridor chat, we 
have noted that it is very tempting to go about trying to prove Conjecture 
15.,'ll by abstract and “easy” arguments, here meaning arguments that do not 
involve any specifics about the geometry or graph structure of To warn 
against such attempts, Deijfen and Hiiggstrom (2006:2) constructed graphs 
on which the two-type Richardson model fails to exhibit the monotonicity 
behavior predicted in Conjecture 15.31 Let us briefly explain the results. 

The dynamics of the two-type Richardson model can of course be defined 
on graphs other than the Z'^ lattice. For a graph Q, write Coex(^) for the 
set of all A > 1 such that there exists a finite initial conhguration (^ 1 ,^ 2 ) 
for which the two-type Richardson model with infection intensities 1 and A 
started from (^ 1 ,^ 2 ) yields infinite coexistence with positive probability. Note 
that, by time-scaling and interchange of the infections, coexistence is possible 
at parameter value A if and only if it is possible at A“^, so no information is 
lost by restricting to A > 1. In Deijfen and Haggstrom (2006:2) examples of 
graphs G are given that demonstrate that, among others, the following kinds 
of coexistence sets Coex(^) are possible: 

(i) Coex(t/) may be an interval (a, 5) with 1 < a < 5. 

(ii) For any positive integer k the set Coex(G) may consist of exactly k 
points. 

(iii) Coex(G) may be countably infinite. 

All these phenomena show that the monotonicity suggested in Conjecture 15.31 
fails for general graphs. However, a reasonable guess is that Conjecture 15.31 is 
true on transitive graphs. Indeed, all counterexamples provided by Deijfen and 
Haggstrom are highly non-symmetric (one might even say ugly) with certain 
parts of the graph being designed specihcally with propagation of type 1 in 
mind, while other parts are meant for type 2. We omit the details. 

6 Unbounded initial configurations 

Let us now go back to the Z'^ setting and describe some results from our 
most recent paper, Deijfen and Haggstrom (2007), concerning the two-type 
model with unbounded initial configurations. Roughly, the model will be 
started from configurations where one of the infections occupies a single site 
in an inhnite “sea” of the other type. The dynamics is as before and also the 
question at issue is the same: can both infection types simultaneously infect 
infinitely many sites? With both types initially occupying infinitely many 
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sites the answer is (apart from in particularly silly cases) obviously yes, so 
we will focus on configurations where type 1 starts with infinitely many sites 
and type 2 with finitely many - for simplicity only one. The question then 
becomes whether type 2 is able to survive. 

To describe the configurations in more detail, write {xi,... ,Xd) for the 
coordinates of a point x E Z'^, and define % = {x \ xi = 0} and C = {x : 
xi < 0 and Xi = 0 for i = 2,..., d}. We will consider the following starting 
configurations. 

liTi) : all points in ?^\{0} are type 1 infected and 

0 is type 2 infected, and . . 

/(£) : all points in T\{0} are type 1 infected and ^ 

0 is type 2 infected. 

Interestingly, it turns out that the set of parameter values for which type 2 is 
able to grow indefinitely is slightly different for these two configuration. First 
note that, as before, we may restrict to the case Ai = 1. Write and q 
for the distribution of the process started from I{'H) and I{C) respectively 
and with type 2 intensity A. The following result, where G 2 denotes the event 
that type 2 grows unboundedly, is proved in Deijfen and Haggstrom (2007). 

Theorem 6.1 For the two-type Richardson model in d > 2 dimensions, we 
have 

(a) P^’q(G 2 ) > 0 i/ and only if X > 1; 

(b) 0 ( 1 ^ 2 ) > 0 i/ and only if X > 1. 

In words, a strictly stronger type 2 infection will be able to survive in both 
configurations, but, when the infections have the same intensity, type 2 can 
survive only in the configuration 

The proof of the if-direction of Theorem lb.ll fa) is based on a lemma stating 
roughly that the speed of a hampered one-type process, living only inside a 
tube which is bounded in all directions except one, is close to the speed of an 
unhampered process when the tube is large. For a two-type process started 
from this lemma can be used to show that, if the strong type 2 infection 

at the origin is successful in the beginning of the time course, it will take 
off along the xi-axis and grow faster than the surrounding type 1 infection 
inside a tube around the xi-axis, thereby escaping eradication. The same 
scenario - that the type 2 infection rushes away along the xi-axis - can, by 
different means, be proved to have positive probability in a process with A = 1 
started from I{C). Infinite growth for type 2 when A < 1 is ruled out by the 
key proposition from Haggstrom and Pemantle (2000) mentioned in Section 
3. Proving that type 2 cannot survive in a process with A = 1 started from 
1(77) is the most tricky part. The idea is basically to divide Z'^ in different 


14 


levels, the l-th level being all sites with xi-coordinate I, and then show that 
the expected number of type 2 infected sites at level I is constant and equal to 
1. It then follows from a certain comparison with a one-type process on each 
level combined with an application of Levy’s 0-1 law that the number of type 
2 infected sites at the Lth level converges almost surely to 0 as Z —)■ oo. 

Finally we mention a question formulated by Itai Benjamin! as well as 
by an anonymous referee of Deijfen and Haggstrom (2007). We have seen 
that, when A = 1, the type 2 infection at the origin can grow unboundedly 
from I{C) but not from It is then natural to ask what happens if we 

interpolate between these two configurations. More precisely, instead of letting 
type 1 occupy only the negative xi-axis (as in I{C)), we let it occupy a cone 
of constant slope around the same axis. The question then is what the critical 
slope is for this cone such that there is a positive probability for type 2 to 
grow unboundedly. That type 2 cannot survive when the cone occupies the 
whole left half-space follows from Theorem l6.ll as this situation is equivalent to 
starting the process from I{l-L). It seems likely, as suggested by Itai Benjamin!, 
that this is actually also the critical case, that is, infinite growth for type 2 
most likely have positive probability for any smaller type 1 cone. This however 
remains to be proved. 
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